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Abstract 



We survey recent results on some one- and two-dimensional patterns generated by random 
permutations of natural numbers. In the first part, we discuss properties of random walks, 
evolving on a one-dimensional regular lattice in discrete time n, whose moves to the right or 
to the left are induced by the rise-and-descent sequence associated with a given random 
permutation. We determine exactly the probability of finding the trajectory of such a 
permutation-generated random walk at site X at time n, obtain the probability measure 
of different excursions and define the asymptotic distribution of the number of "U-turns" 
of the trajectories - permutation "peaks" and "through". In the second part, we focus 
' on some statistical properties of surfaces obtained by randomly placing natural numbers 

1, 2, 3, . . . , L on sites of a Id or 2d square lattices containing L sites. We calculate the 
distribution function of the number of local "peaks" - sites the number at which is larger 
than the numbers appearing at nearest-neighboring sites - and discuss some surprising 
collective behavior emerging in this model. 



1 INTRODUCTION 



Properties of random or patterns-avoiding permutations of series of natural numbers have 
been analyzed by mathematicians for many years. Studies of several problems emerging within 
this context, such as, e.g., the celebrated Ulam's longest increasing subsequence problem (see, 
5— i e.g., Refs. ^GlEl and references therein), provided an entry to a rich and diverse circle of 
mathematical ideas 4 , and were also found relevant to certain physical processes, including 
random surface growth [SHEHS] or 2D quantum gravity (see, e.g., Ref. 0). 

In this paper we focus on random permutations from a different viewpoint addressing the 
following question: what are statistical properties of patterns, e.g., random walks or surfaces, 
when random permutations are used as their generator? We note that such a generator is 
different of those conventionally used, since here a finite amount of numbers is being shuffled 
and moreover, neither of any two numbers in each permutation may be equal to each other; 
this incurs, of course, some correlations in the produced sequences of random numbers. 
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In the first part of this paper, we consider a simple model of a permutation-generated 
random walk (PGRW), first proposed and solved in Ref. UJ. In this model, random walk 
evolves in discrete time on a one-dimensional lattice of integers, and the moves of the walker 
to the right or to the left are prescribed by the rise-and-descent sequence characterizing each 
given permutation ir — {tti, ^2,^3, 7T n +i}. In a standard notation, the "rises" (the 

"descents") of the permutation n are such values of I for which 717 < 717+1 {^1 > Tj+i). 

We determine exactly several characteristic properties of such a random walk, including the 
probability V n (X) of finding the end-point X n of the PGRW trajectory at site X, its moments, 
the probability measure of different excursions and the asymptotic distribution of the number 
of the "U-turns" of the trajectories, which, in the permutation language, corresponds to the 
number of "peaks" and "through" in a given permutations. 

In the second part, we focus on some statistical properties of surfaces created by randomly 
distributing numbers 1,2,3, ... ,L on sites of one- or two-dimensional lattices containing L 
sites. A number appearing at the site j determines the local height of the surface. Denoting 
as local surface "peaks" such sites the number at which exceeds the numbers appearing at 
neighboring sites, we aim to calculate the probability P(M, L) of having M peaks on a lattice 
containing L sites. 

In one dimension this can be done exactly ^H]- We are also able to calculate the "correlation 
function" p(l) defining the conditional probability that two peaks are separated by the interval 
/ under the condition that the interval / does not contain other peaks. In 2D, determining 
exactly first three cumulants of P(M, L), we define its asymptotic form using expansion in the 
Edgeworth series and show that it converges to a normal distribution as L — > 00 10 . For 
2D model, we will also discuss some surprising cooperative behavior of peaks. 

2 PERMUTATION-GENERATED RANDOM WALK 

Suppose there are two players - A and B going to play the following mindless card game. 
They first agree on the value of each card linearly ordering them, e.g., by putting suits in the 
bridge- bidding order: Clubs < Diamonds < Hearts < Spades, and adopting a convention that 
in a series of cards with the same suit "two" is less than "three", "three" is less than "four" 
and etc, while the ace has the largest value. So, in such a way, a deck of cards is labeled 
1,2,3,..., 52. Then, the deck is shuffled, the upper card is turned its face up, and our players 
start the game: the second card is turned face up; if its value is higher than the value of the 
first card, the player A receives some unit of money from the player B\ if, on contrary, its value 
is less than the value of the first card, player B receives a unit of money from the player A. 
At the next step, the third card is turned its face up and its value is compared against the 
value of the second; if, again, its value is higher than the value of the second card, the player 
A receives money from the player B, otherwise, the player A pays the player B. The process 
continues until the deck is over. One is curious, as usual, about the winner and the amount of 
his gain. 

Let us now look on such a random process more formally. Let 7r = [711, it2, 1T3, ■ . ■ , ^n+i] 
denote a random permutation of [n + 1]. We rewrite it next in two- line notation as 



and suppose that this table assigns a discrete "time" variable I (I = 1, 2, .3, n + 1, upper line 
in the table) to each permutation in the second line. We call, in a standard notation, as a "rise" 
(or as a "descent") of the permutation tt, such values of I for which 7T; < -ki+i (7r2 > tti+i). 




(2.1) 
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Consider now a Id lattice, place a walker at the origin at time moment I = 0, and let it 
move according to the following rules: 

-at time moment I = 1 the walker is moved one step to the right if tti < 7r 2 , i.e., I = 1 is a rise, 
or to the left if 7i"i > ir 2 , i.e I = 1 is a descent. 

-at I = 2 the walker is moved to the right (left) if TT2 < Tt3 (tV2 > ^3, resp.) and etc. 

Repeated I times, this results in a random trajectory x\ n \ (I = 1, 2, . . . , n), such that 

X i n) =i>2 s k, Sfe = sgn(7r fc+ i -7r fc ), sgn(x) = \ + ] \f >0 ' (2.2) 
f— ; —1 otherwise. 

Evidently, .Xj is just an amount of money the player A won (or lost, if x[ n ^ < 0) up to 
the moment when the l-th card is opened. Below we answer on a set of questions on various 

(n) 

statistical properties of random variable X t . 

2.1 Probability Distribution of The End-Point of Trajectory xj n '. 

Let V n (X) denote the probability that the walker is at site X at time moment / = n. Let N\ 
(A/j.) be the number of "rises" ("descents") in a given permutation ir. Evidently, the end-point 
X n of the walker's trajectory is X n = A/f — A/j_. Since A/f +A/j. = n, we have that X n = 2A/j — n 
and hence, X„ is fixed by the number of rises in this permutation. 

Number of permutations of with exactly A/f rises is given by the Eulerian number |12j : 

n ;, , )-|V l r(»;>, + .-r,-, (r)-^^- m 

Consequently, V n {X) is given explicitly by [5j 

v m _ [1 + (-1) X+ "1 /n+l\ , . 

- 2(n + l)! \ / ' (2 ' 4) 

Several useful integral representations of the distribution function V n (X) and of the corre- 
sponding lattice Green function were also derived |3J . Expression in Eq. (|2.17|l may be cast into 
the following form: 

„. w _ n +(-')"•] r «*n m(Xk)dk . (2 . 6) 



Note that it looks almost like usually encountered forms for Polya walks JSIi just the upper 
limit of integration is infinity but not "7r"; as a matter of fact, an integral representation with 
the integration extending over the first Brillouin zone only obeys 

and has a very different structure compared to that of Polya walks. So does the lattice Green 
function Q(X,z), associated with the result in Ea. l|2.t)[) : 

q{x, z) =± v n (x) Z " = J- r ^p^m dk . (2 . 7) 

nz J sin (k — zsm{k)) 
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To establish the asymptotic form of V n (X) we go to the limit z — * 1 , and, inverting the 
expression in Ea. (|2.7|) in this limit, we find that V n (X) converges to a normal distribution: 

Consequently, in this limit the correlations in the generator of the walk - random permutations 
of [n+1], appear to be marginally important; that is, they do not result in anomalous diffusion, 
but merely affect the "diffusion coefficient" making it three times smaller than the diffusion 
coefficient of the standard ID Polya walk. 



2.2 Correlations in PGRW Trajectories. 

So, we have a symmetric random walk, which makes a move of unit length at each moment of 
time with probability 1, but nonetheless looses somehow two thirds of the diffusion coefficient. 
To elucidate this puzzling question, we address now an "inverse" problem: given the distri- 
bution V n (X), we aim to determine two- and four-point correlations in the rise-and-descent 
sequences and consequently, in the PGRW trajectories. 
For the Fourier-transformed V n {X) we find 



V n {k)= V exp(zkX)T n (X)= [ > ; Li- w - 1 (cxp(-2ifc)), 



(2.9) 



Li_„_i(y) being a poly-logarithm function. Ea. l|2.9|l yields 

Xfy = i(n+2-25 n>0 ) and (xfy = ^((5n+8)(n+2)-165 n ,o-245 nil +85 n , 2 ), (2.10) 
where &y is the Kronecker delta, such that their generating functions are 



n=0 3 t 1 - Z ) 



n— ^ ' 

On the other hand, the second and the fourth moments of the walker's displacement obey 

-I 2 X 



(Xl) = ((A/- T -A/i) 2 ) = 



E 

.1=1 



Si 



n—1 n 

(2) 



E <£ 

jl=l 32=71 + 1 



<A^) = n + n( " 2 lj fl( n - 2) + (8(n - l)fl(n - 2) + 12 (n(n - 3) + 2) 0(n - 3))c (2) (l) + 

+ 24 V (n-2-m)e^(m)0(n-4), 8{x) = j +1 lf 1 > °' (2.12) 
z — ' V / otherwise. 

m— 1 ^ 

where C^ 2 ^(m) = cj 2 ^- 1 = (sj 1 Sj m+1 ) is the two-point correlation function of the rise-and- 
descent sequence, while C^(m) = ( S;/l s ii+i s ji+m+i s ji+m+2) defines the four-point correlation 
function. Note that fc-point correlations with k - odd vanish [S] and thus C^\m) is an only 
non- vanishing form of four-point correlations. 
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Multiplying both sides of Eas. (|2.12|l by z n and performing summation, we find that 
V ) ^ V J 2z w 2 3' 

m— 1 

- 16z 2 - 7z - 3 



24z 3 v ' 72z 2 (l-z) 



which imply 



m— 1 



c(2) m= n 1/3 m ; 2 c(4) ^= w tf f m >J (2 - w) 

10 m > 2 11/9 if m > 2 



Equations (|2.14() signify that for m > 2 the two-point correlations C^ 2 ^(m = \j% — in 
the rise-and-descent sequences decouple into the product < Sj 1 >< Sj 2 > and hence, vanish. 
Consequently, two-point correlations extend to the nearest neighbors only. The four-point 
correlations decouple C^{m) =C' 2 '(l)C' 2 '(l)form>2. 

Equations (|2.14|) enable us to calculate explicitly the probabilities of different configurations 
of two, three and four rises and descents, which are listed below 

J 1/6 m = l , . J 1/3 m = l 

I 1/4 m > 2 I 1/4 m > 2 

Hence, configurations with two neighboring rises or two neighboring descents have a lower 
probability than mixed rise-descent or descent-rise sequences; in other words, rises "repel" 
each other and "attract" descents. Further on, we have: 

, . (l/4 m=l , . fl/24 m = l , , f 1/120 m = l 

PTu( m ) = < 1 /10 . PTT,t("») = S 1 . _ Pn,n( m ) = < , ,„ R . 

I 1/ 12 m>2 I ^ m>2 [l/'JO m>2 

/3/40 m=l fl/30 ? 7i=l 

11/18 m>2 11/18 ?ti > 2 



. . J 2/15 m = \ , , J 1/20 m = l . , J 11/120 m = l 

^.uH = | 1/9 PTT,uM = | 1/36 m ^ 2 WT,u(m) = | 1/g ^ 

Analyzing these results, we notice that the probabilities of different rise-and-descent sequences 
depend not only on the number of rises or descents, but also on their order within the sequence. 
Correspondingly, reconstructing the PGRW trajectories Xj with n — 4 (see Fig.l) we notice 
that the PGRW transition probabilities depend not only on the number of steps to the right 
or to the left, which the walker has already made, but also on their order. In other words, the 
PGRW represents a genuine non-Markovian process with a memory. 

Finally, for fc-point correlation functions = (sj 1 Sj 2 sj 3 ■ ■ ■ sj k ) of the rise-and-descent 

sequences we find the following behavior [H] : C^j j k factorizes automatically into a product of 
the corresponding correlation functions of the consecutive subsequences, in which all jk differ 
by unity, as soon as any of the distances jk+i — jk exceeds unity. On the other hand, the 
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-5-4-3-2-1 1 2 3 4 



X 

Fig. 1. A set of all possible PGRW trajectories Xf for n = 4. Numbers above the solid arcs 
with arrows indicate the corresponding transition probabilities. Dashed-lines with arrows connect the 
trajectories for different values of I. 



correlation function = - +k of a consecutive sequence of arbitrary order k can be 

obtained exactly. One has that obey the following recursion 



C W = ^ii)^ C {W-P) ) C(°) = l, (2.15) 
p =o + !)■ 

which implies that are related to the tangent numbers |15j . and are given explicitly by 

(_l)fc 2 *+ 2 (2*+ 2 - 1) , s 

Ci ) = ( fe + 2)! LB ™> ^ 

where B/~ are the Bernoulli numbers. Note that since the Bernoulli numbers equal zero for k 
odd, correlation functions of odd order vanish, i.e., C( 2fc +t) = 0. Ea. l|2.16(l also implies that for 
k 3> 1, the fc-point correlation functions (k is even) decay as ~ 2(— \) k / 2 [2/i:) k+2 . 

2.3 Probability Distribution of x( n) for 1 < n. 

(n) 

To determine the structure of excursions Xf of the PGRW we adapt a method proposed by 
Hammersley ^3] m his analysis of the expected length of the longest increasing subsequence. 
The basic idea behind this approach is to build recursively an auxiliary Markovian stochastic 
process Yi, which is distributed exactly as Xj. 

At each time step I, we define a real valued random variable xi+i, uniformly distributed 
in [0,1]. Further on, we consider a random walk on a one-dimensional lattice of integers 
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whose trajectory Yj is constructed according to the following step-by-step process: at each 
time moment I a point-like particle is created at position xi+%. If xi + \ > Xi, a walker is moved 
one step to the right; otherwise, it is moved one step to the left. The trajectory Yj is then 
given by Y; = J2k=i s E n ( x k+i — %k), where sgn(x) is defined in Ea. (|2.2l) . 

We note that the joint process (xj+i, Yj), and therefore Y;, are Markovian since they depend 
only on [xi, Yj_i). Note also that Y; is the sum of correlated random variables; hence, one has 
to be cautious when applying central limit theorems. A central limit theorem indeed holds for 
the Markovian process Y;, but the summation rule for the variance is not valid. 

Two following results have been proven in Ref. [5| concerning the relation between this 
recursively constructed Markovian process Y; and the PGRW: 

(a) the probability P(YJ = X) that the trajectory Y; of the auxiliary process appears at site 
X at time moment I is equal to the probability Vi(X) — P(Xj = X) that the end-point Xi 
of the PGRW trajectory X^ generated by a given permutation of [I + 1] appears at site X. 

(b) the probability P(XP^ — X) that the PGRW trajectory X^ at intermediate time moment 
I = 1, 2, 3, . . . , n appears at site X is equal to the probability P(Yj = X) that the trajectory of 
the auxiliary process Y/ appears at time moment / at site X. In other words, one has 

P(Y = X) = P(X,W = X) = V t (X) = [ ^I0P (^i 1 ) ■ (2-17) 

in) 

Note that distribution of any intermediate point A z of the PGRW trajectory generated by 
permutations of a sequence of length n + 1 depends on I but is independent of n. 



2.4 Probability Measure of PGRW Trajectories. 

(n) 

The equivalence of the processes Y\ and X\ allows to determine the probability measure of 
any given trajectory. We note that in the permutation language, this problem amounts to 
the calculation of the number of permutations with a given rise-and-descent sequence and has 
been already discussed using an elaborated combinatorial approach in Refs. [18116117] . A novel 
solution of this problem has been proposed in Ref. 9 , which expressed the probability measure 
of any given PGRW trajectory (or of some part of it) as a chain of iterated integrals. 
Consider a given rise-and-descent sequence a{k) of length k of the form: 

«(*> = ( 1 2 3 - M. 

V a i o 2 a 3 ... a k J 

where ai can take either of two symbolic values — j or J,. Consequently, the first line in the 
table is the running index I which indicates position along the permutation, while the second 
line shows what we have at this position - a rise or a descent. Assign next to each symbol at 
position / an integral operator; h{\) for a rise (f) and /;(!) for a descent {{): 

/i(T)= / dxi and Ii{l)= f ' 1 d Xl . (2.18) 

Define next a characteristic polynomial Q(x,a(k)) as an "ordered" product 9 



k 

Q(x,a(k)) = :JJ: h{ai) ■ 1, 
i=i 



(2.19) 
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where a; = {T,i} for I = l,...,k. The probability measure P(a(k)) of this given rise-and- 
descent sequence a(k) in the ensemble of all equally likely permutations is then given by 



P(a(k)) 



Q(x, a(k))dx. 



(2.20) 



In may be expedient to illustrate this formal consideration on a particular example. Let a 
given rise-and-descent sequence be of the form {1\T>J.)t>t}- For this sequence, the character- 
istic polynomial Q(x, ct(5)) and the probability of this particular configuration obey 



Q(a? ) o(5))=7i(t)/2(T)ia(i)/4(T)iB(T)-l= / dxi / dx 2 / dx 3 / dx 4 / dx 5 • 1 



3 

40 



x 
8 



x 
12 



x 
24 



120 



f 19 

and P(o(5)) = y Q(x,a(5)) dx = —. 



(2.21) 



Another way to look on the problem of calculation of the measure of a given PGRF tra- 
jectory is to use the results of Niven determining the number of permutations with a given 
rise-and-descent sequence |17j . Following Niven, consider a fixed up-and-down arrow sequence 
a(k) of length k and denote by l\, I2, ■ ■ ■ , l r the positions of downarrows (descents), r being the 
total number of downarrows along the sequence. A question now is to calculate the number 
N(Xi) of permutations generating a given up-and-down sequence (or, in our language, a 

given trajectory A^™" 1 ). Combinatorial arguments show that N(Xt ) equals the determinant 
of a matrix of order r + 1 whose elements ciij (where i stands for the row, while j - for the 
column) are binomial coefficients ( ; 4 1 ) , where Iq — and l r +i = k + 1 17 . Consequently, an 
alternative expression for the probability P(a{k)) may be written down as 



(2.22) 



One may readily verify that both Eg. (|2. 20(1 and Ea. l2.22() reproduce our earlier results de- 
termining the probabilities of different four step trajectories. Note also that the probability 
measure defined by Ea. H2.20[) or Ea. (|2.22|l is not homogeneous, contrary to the measure of the 
standard Polya walk. 
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2.5 Distribution of The Number of "U-Turns" of PGRW Trajectories. 

In this subsection we study an important measure of how scrambled the PGRW trajectories 
are. This measure is the number Af of the "U-turns" of an rt-step PGRW trajectory, i.e. the 
number of times when the walker changes the direction of its motion. 

In the permutation language, each turn to the left (right), when the walker making a 
jump to the right (left) at time moment I jumps to the left (right) at the next time moment 
I + 1 corresponds to a peak f| (a through JT) of a given permutation tt, i.e. a sequence 
iri < 7T; + i > 7T; + 2 (j^i > ^1+1 < 7r i+2)- Consequently, the distribution function V(Af,n) of the 
number of the "U-turns" of the PGRW trajectory, (i.e. the probability that an n-step PGRW 
trajectory has exactly Af turns), is just the distribution of the sum of peaks and through. 
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Let Af p and Aft denote the number of peaks and through in a given random permutation 
[n + 1] , respectively. These realization-dependent numbers can be written down as 



ra-l 



n-1 



Af P =^(l + Sl )(l- %1 ), M = ^(l-sO(l + sw), (2-23) 
1=1 i=i 

such that total number Af of U-turns of the PGRW trajectory is given by 

_^ n— 1 

Af = Af p +Af t = - ^ (1 - 



(2.24) 



One finds then that the characteristic function Z n (k) of A/" is a polynomial in tanh(z/c/2) 9 : 



Z n (k) = ^ exp (ikAf) j = exp 
[n/2] / 



exp 



ifc(n — 1) 



— 2Z + 1\ (£j=i 



Z=l 



1 + e 1 



\ Y"! ■ i m> ; / mi !m2 ! ... to/! 

,lm 1 +2m 2 + ...+im ! =i V ^3=1 J 7 j=l 

where C^ 2 -?) obey Eq.(2.16). The generating function of Z n (k) is given by 

4 2 

Z„(fc) z" =-- 

r 1 _i_ 

n=2 



Z(fc,z) = ^Z„(fc) 



(l + e< fc ) z l + e lfc 



z + 



(l + e lfc ) z 
4 

(l + e ife ) 2 z 



3=0 



1-e 



2 



1 + e 



coth ( 1 - 



' (l — e 2ik ) z 2 \ 
^•) 1/2 |)-1 



1 + e 



Turning next to the limit z — > 1 _ and inverting Z(fc, z) with respect to fc and z, we find that in 
the asymptotic limit n — > oo, the distribution function V(Af, n) of the number of the "U-turns" 
in the PGRW trajectory converges to a normal distribution: 



3/ 5 \i/2 



/ 



exp 



V(Af,n)^U^) 
with mean value 2n/3 and variance a 2 = 8n/45 



45(AT- |n 
16n 



(2.25) 



2.6 Diffusion Limit 

Consider finally a continuous space and time version of the PGRW in the diffusion limit. In 
Ref. UJ, the following approach has been developed: 

Define first the polynomial: V"'(x,Y) — J2yi=y Q( x > a (0); where the polynomial Q(x,a(l)) 
has been determined in Eq. (|2.19|l and the sum extends over all Z-steps trajectories Y\ starting at 
zero and ending at the fixed point Y. Note next that one has Vi(Y) — (V^{x, Y)} , , . Now, 
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for the polynomials V®(x,Y) one obtains, by counting all possible trajectories Y. starting 
from zero and ending at the fixed point Y , the following "evolution" equation: 

V^ l+1 \x,Y) = J,(T) • V^(x,Y- 1) + 1,(1) • V^(x,Y + 1). (2.26) 

(n) 

Taking next advantage of the established equivalence between the processes YJ and Y ; , we 
can rewrite the last equation, upon averaging it over the distribution of variables {x{\, as 

V l+l {Y) = - 1) + l^f l{ Y + 1), (2.27) 

which represents the desired evolution equation for Vi(Y) in discrete space and time. 

We turn next to the diffusion limit. Introducing y = aY and t = rl variables, where a and 
t define characteristic space and time scales, we turn to the limit a, t — > 0, supposing that the 
ratio a 2 /r remains fixed and determines the diffusion coefficient Dq = a 2 /2r. In this limit, 
Ea. H2.27|) becomes 



Note that the resulting continuous space and time equation is of the Fokker-Planck type; it 
has a constant diffusion coefficient and a negative drift term which, similarly to the Ornstein- 
Uhlenbeck process, grows linearly with y, but the amplitude of the drift decays in proportion 
to the first inverse power of time, which signifies that the process Y, eventually delocalizes. 
The Green's function solution of Ea. (|2.28|) . remarkably, is a normal distribution 



which is consistent with the large-n limit of the discrete process derived in Ea. H2.8|) . 



3 RANDOM SURFACES GENERATED BY RANDOM PERMUTATIONS. 

Consider a one-dimensional lattice containing L sites on which we distribute at random num- 
bers drawn from the set 1, 2, 3, . . . , L, (see Fig[2J). Suppose that a number appearing at the site 
j determines the local height of the surface. Further on, we denote a local "peak" as a site j the 
number at which exceeds the numbers appearing at two neighboring sites. Generalization to 
two-dimensional square lattice with L = m x m, (m is an integer) sites is straightforward (see 
FigEJ : the only difference here is that we call as local surface peaks such sites j the numbers 
at which are greater than numbers appearing at four adjacent sites. 

For these models, our goal is to evaluate the probability P(M, L) that the surface created 
in such a way has exactly M peaks on a lattice containing L sites |10| . In one dimension this 
can be done exactly and provides also a distribution function of the number of right U-turns 
of the PGRW trajectories. Here we are also able to calculate the "correlation function" p(l) 
defining the conditional probability that two peaks are separated by the interval / under the 
condition that the interval I does not contain other peaks. In 2D, determining exactly first 
three cumulants of P{M, L) we define its asymptotic form using expansion in the Edgeworth 
series JT] and show that it converges to a normal distribution as L — > oo ^Hj- For 2D model, 
we will also discuss some surprising cooperative behavior of peaks. 
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Fig. 2. (a) One-dimensional (L — 6) and (b) two-dimensional square (L = 4 x 4) lattices with 
periodic boundary conditions on sites of which we distribute numbers drawn from the set 1, 2, 3, . . . , L. 
These numbers determine the local heights of the surface. The sites with numbers bigger than those 
appearing on the neighboring sites are refereed to as the "peaks". 



3.1 Probability Distribution P(M,L) in ID Case 

In one dimension, the probability P{M, L) of having M surface peaks on a lattice containing 
L sites can be calculated exactly through its generating function [10) : 

oo 

W(s,L)=L\ s M+1 P(M,L), (3.1) 

M=0 

where 

Inverting ()3.1[l . we find the following exact expression for P(M, L): 

oL+2 M / L+l \ ' m i + 1 

One finds then that in the limit L — * oo, the probability P(M,L) converges to a normal 
distribution: 

3 HT f 45(M-±L) 2 



p ^>~2V^ exp r 4l j i (3 - 4) 

with mean L/3 and variance cr 2 = 2L/45. Note that this result is consistent with the distribu- 
tion of the number of U-turns of the PGRW trajectories, Eq.(2.25). 

3.2 Conditional Probability p(l) of Two Peaks Separated by Distance I 

We aim now at evaluating the conditional probability p(l) of having two peaks separated by 
a distance I, under the condition that there are no peaks (i.e. sequences 1 1) on the interval 
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between these peaks. Using the approach put forth in Section 2.4, we have that this probability 
obeys 



P(l) 



J dx^2Q(x,a), 



(3.5) 



where the sum is taken over all possible peak-avoiding rise-and-descent patterns of length I 
between of two peaks, while Q(x,a) denote the Q-polynomials corresponding to each given 
configuration. 

There are several possible peak-avoiding rise-and-descent sequences contributing to such a 
probability. These sequences are depicted in Fig. 3. 
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Fig. 3. Rise-and-descent patterns contributing to the conditional probability of having two closest 
peaks at distance I apart from each other. Configuration (a) has only one descent between two peaks. 
Configuration (b) has two descents following the first peak, i.e. a "through" at xi, and (c) presents 
a generalization of (b) over configurations having m descents, (m = 1,2,...,/), after the first peak 
which are followed by I — in rises, i.e. a "through" at x m -i- 



The Q-polynomial associated with the general configuration (c) in Fig. 3 is given by 

i-\ rX ,-Y i-X! rXm-2 rl rl rxi-i 

Q c ( x ) = / dX / dY / dxi / dx 2 ■ • • / dx m -\ / dx m . . . dxi-i / dxi 

Jx JO JO Jo Jo ix m -i Jxi-2 Jo 

(3.6) 

Calculating this integral recursively and summing over m, m — 1,2, ... ,1, we find that the 
generating function of the probability p{l) obeys 

1=2 v 

2 12 1 4 

15 9 35 45 567 ^ ' 
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Inverting the latter expression, we find that p(l) is determined explicitly by 



p(0 



2 i+3 



2 \(l + iy. (Z + 2)! (1 + 3)! 



2 i+i 2 l + 2 
- 2- 



= 2 



t (Z -l)(f + 2) 
(i + 3)! : 



(3.8) 



Note also that p(l) in ea. (|3.8(l coincides with the distribution function of the distance between 
two "weak" bonds obtained by Derrida and Gardner in their analysis of the number of 
metastable states in a one-dimensional Ising spin glass at zero temperature. 

3.3 Probability Distribution P(M,L) in 2D Case 

We turn now to calculation of the probability P(M, L) of having M peaks on a two-dimensional 
square lattice containing m x m = L sites. We linearly order lattice sites by index j, j = 
1, 2, . . . , L, in the same way as an electron beam highlights the TV screen. Note that in such a 
representation a site j is a peak if and only if 7Tj is simultaneously larger than TTj—i, ffj+i, ftj-m 
and Wj+ m - 



(a) 



(b) 
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l(x) - normal distribution 
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-4-3-2-101234 
X 

Fig. 4. (a) Results of numerical simulation of p(x, L) for L — 100, 256: (a) comparison of p(x, L) with 
the Gaussian function g(x); (b) comparison of the function p(x, L)/ g(x) with 1 + /(x)/vL. 



Using the approach devised in Section 2.4, we find the first cumulant of P(M, L): 



„■ Jo Jo u 



kO(xj - Xj-\)6{xj - Xj+i)6{xj - x j+m )9(xj - Xj- m ) = -L (3.9) 



where the summation over "j" and "k" extends over all lattice sites, i.e. j,k = 1,2,3,..., L. 
Further on, we find that the second cumulant is given by 



2D 



EE 

j » 



Y^dxkOixj - Xj..i)6(xj - x j+ i)0(xj - x j+m )0(xj - x^ m ) x 

(3.10) 



1 13 

x 9(xi — Xi-i)9(xi — Xi + i)9(xi — Xi +m )9(xi — Xi- m ) — — 225"^ 
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Calculation of the third cumulant is more involved JUj and yields 

4° = ^L. (3.11) 

As a matter of fact, we can show that cumulants of any order grow in proportion to the 
first power of L. Introducing next the normalized deviation, x = (M — k\ d )/ n\ D , we seek 

the normalized probability distribution p(x, L) = P Lj expanding it in the 

Edgeworth series (cumulant expansion) We find then that p(x,L) obeys 

p(x, L) ~ g(x) (l + J=/(x) + o L+=\ \ , (3.12) 



(3.13) 



where g{x) is the Gaussian function g(x) = exp(— x 2 /2)/2ir and f(x) is given by 

rr K 3° !/ <\ o s 512 /225\ 3/2 l, , „ s 
^=^7^6^ " 3X) = 32175 ("13 J 6^ - 3 ^- 



Note that f(x) is independent of L and correction terms in Ea. H3.12|l decay faster than 1/y/L. 

To check the asymptotical result in Ea. (|3.12|l . we have performed numerical simulations 
for the discrete 2D permutation-generated model with periodic boundary conditions and have 
computed the distribution function p(x, L) numerically. In Fig. 4a we present numerical data 
for p{x, L) and compare it against the Gaussian function g{x) for system sizes L = 100 and 256. 
Furthermore, in Fig. 4b we plot the ratio p(x, L)/g(x) as the function of x. One notices that 
the deviation of the numerically computed function p(x, L) from the Gaussian function g{x) 
is indeed very small. Moreover, the difference between the normalized probability distribution 
function p{x 1 L) and the Gaussian function g(x) is smaller the larger L is. 



3.4 Gas, Liquid and Solid of Peaks 

In this last subsection we discuss a surprising collective behavior of peaks. Consider a somewhat 
exotic "statistical physics" model with the partition function: 

all permutations 

where z is a fugacity and M determines the number of peaks for a given distribution of natural 
numbers 1, 2, 3, . . . , L on a lattice containing L sites. In one- and two-dimensions 

3 

M = ^ °( 7r 3 - Kj-i WiTTj - TT j+ i)9(nj - nj +m )0{iTj - nj- m ) (3.15) 

3 

respectively. A striking feature of this model in two-dimensions, as communicated to us by 
B.Derrida who analyzed a similar partition function numerically, is that it undergoes a phase 
transition at z « exp(2.2) |20|! 

Below we present some arguments why such a transition may indeed take place in two- 
dimensions: Note first that peaks, by definition, can not simultaneously occupy nearest- 
neighboring lattice sites. It means that if we view the peaks as some fictitious particles, they 
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Fig. 5. Probabilities of having a single peak and a pair of isolated peaks some distance apart from 
each other on one- and two-dimensional lattices. Filled circles denote peaks - sites with numbers ex- 
ceeding numbers appearing on nearest-neighboring sites. First column presents probabilities of several 
configurations in Id, while the second one describes the probabilities of a few possible configurations 
appearing in two dimensions. 



are not point particles occupying just a single site but rather "hard-squares". In principle, 
this is already enough to produce a phase transition, but at somewhat higher values of z |21j , 
since here a cost of each particle is higher requiring a particular arrangement of numbers on 
five lattice sites. 

Still striking, it appears that peaks experience short-range "attractive" interactions. Note 
that a probability of having a single isolated peak is equal to 1/5. The probability of having 
two peaks separated by distance exceeding I = 2 is 1/25 = 1/5 2 , which means that at such 
distance the peak do not feel each other. On the other hand, the probability of having two 
peaks at next-nearest-neighboring sites or in corners of each lattice cell, is greater than 1/25, 
which means that peaks effectively "attract" each other! In Fig. 4 we summarize the results for 
the probabilities of several simple configurations of peaks. 

This " attractive" force between peaks can be thought off as a sort of a depletion force acting 
between colloidal particles: by definition, each peak is a site with a number which is greater 
than numbers on four neighboring sites. When two peaks are next-nearest- neighboring sites, 
they have a "common" site with a number which is restricted to be less than the least of the 
numbers on the peak sites, i.e. may attain less values than numbers on the sites around an 
isolated peak. This effect is even more pronounced for peaks occupying sites in the corners of 
a lattice cell - here these two peaks have two " common" numbers and each of them has to be 
less than the least of them. 

Note finally that for clusters containing several peaks, the presence of such common sites 
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makes the interactions between peaks non-additive and dependent on the specific features of 
the cluster's geometry. Thus it is not evident at all whether the phase transition is in the Ising 
universality class. 

4 ACKNOWLEDGMENTS 

The authors wish to thank B.Derrida, P.Hanggi, P.Krapivsky and Z.Racz for helpful discussions 
and remarks. The work is partially supported by the grant ACTNIM-2004-243 "Nouvelles 
Interfaces des Mathematiques" (France). 

References 

[1] Baik, J., Deift, P. & Johansson, K. 1999 J. Amer. Math. Soc. 12 1119 

[2] Deutsch, E., Hildebrandt, A.J. & Wilf, H.S. 2003 Elec. J. Combinatorics 9(2) R12 

[3] Odlyzko, A.M. & Rains, E.M. 2000 Contemp. Math. 251 439 

[4] Aldous, D. & Diaconis, P. 1999 Bull. Amer. Math. Soc. 36 413 

[5] Seppalainen, T. 1996 Elec. J. Prob. 1 1 

[6] Prahofer, M. & Spohn, H. 2000 Phys. Rev. Lett. 84 4882 

[7] Johansson, K. 2000 Comm. Math. Phys. 209 437 

[8] Majumdar, S.N. & Nechaev, S. 2004 Phys. Rev. E 69 011103 

[9] Oshanin, G. & Voituriez, R. 2004 J. Phys. A: Math. Gen. 37 6221 
[10] Hivert, F., Nechaev, S., Oshanin, G., & Vasilyev, O. On the distribution of surface extrema in several 

one- and two-dimensional random landscapes J. Stat. Phys., to appear; cond- mat/0509584 
[11] Cramer, H. 1957 Mathematical Methods of Statistics (Princeton University Press: Princeton) 
[12] Graham, R.L., Knuth, D.E., &: Patashnik, O. 1994 Eulerian Numbers, in: Concrete Mathematics: A 

Foundation for Computer Science, (Addison-Wesley, Reading MA) 
[13] Hughes, B.D. 1995 Random Walks and Random Environments Vol.1 (Oxford, Clarendon) 
[14] Hammersley, J.M. 1972 Proc. 6th Berkeley Symp. Math. Statist, and Probability 1 345 (University of 

California Press) 
[15] see, e.g., Foata, D. & Strehl, V., 1974 Math. Z. 137 257 
[16] Carlitz, L. 1973 Math. Nachr. 58 31 
[17] Niven, I. 1968 Nieuw Atchief voor Wiskunde 16 116 

[18] Macmahon, P.A. 1915 Combinatorial Analysis (Cambridge University Press) 
[19] Derrida, B. & Gardner, E. 1986 J. Phys. (Paris) 47 959 
[20] Derrida, B. unpublished 

[21] Baxter, R.J. 1982 Exactly Solved Models in Statistical Mechanics (Academic Press, New York) 



